Field— induced structural aging in glasses at ultra low temperatures 
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In non-equilibrium experiments on the glasses Mylar and BK7, we measured the excess dielectric 
response after the temporary application of a strong electric bias field at mK-temperatures. A model 
recently developed describes the observed long time decays qualitatively for Mylar Q], but fails for 
BK7. In contrast, our results on both samples can be described by including an additional mechanism 
to the mentioned model with temperature independent decay times of the excess dielectric response. 
As the origin of this novel process beyond the "tunneling model" we suggest bias field induced 
structural rearrangements of "tunneling states" that decay by quantum mechanical tunneling. 

PACS numbers: 61.43.Fs, 77.22.-d, 66.35.+a, 05.70.Ln 



Most low temperature properties of glasses are well de- 
scribed by the phenomenological "tunneling model" 0, 
E3- [3- E] ■ It starts from the assumption that the poten- 
tial minima of groups of atoms are not " well defined" , but 
can be described as particles moving in " double well" po- 
tentials. At low enough temperatures, only the ground 
states of the two wells are occupied, and the "tunneling 
states" (TSs) move by quantum mechanical tunneling be- 
tween them. The energy splitting of these two levels is 
given by E = \J Aq + A 2 , where the asymmetry, A, is 
the energy difference between the two local minima. The 
tunneling splitting can be estimated using the WKB- 
approximation as Ao ~ ?if2exp(— A), with the attempt 
frequency, 0. The tunnel parameter A = d\/2mV /2h 
contains the effective mass, m, the potential barrier 
height, V, and the distance between the two wells, d. 
Recent experiments revealed deviations from the tunnel- 
ing model and show that interactions between TSs gain 
importance with decreasing temperature 0, 0, ■ Our 
new results will be discussed within a picture that goes 
beyond "double well" potentials. 

Long time non-equilibrium dynamics of TSs in the 
polyester glass Mylar after the temporary application of 
a large electric bias field were previously presented 0. 
These data were analyzed using the "dipole gap" the- 
ory H by extending its predictions to the case of a tem- 
porary applied bias field. A central assumption of the 
"dipole gap" theory is that mutually strongly coupled 
TSs do not contribute to the dielectric response in ther- 
mal equilibrium. However, out of equilibrium these TSs 
contribute to the dielectric constant as if isolated § . The 
findings in Q were explained in terms of strongly coupled 
pairs of TSs that are driven out of equilibrium during 
a field sweep, not only by relaxation, but also by non- 
adiabatic driving of one member of two strongly coupled 
TS-pair by the applied field. Our new results indicate 
the existence of an additional process leading to non- 
equilibrium pairs of TSs, while the sum of all processes 
remains constant. This novel process yields an excess di- 
electric response with a decay time independent of the 
bias field history and temperature. For the polyester 



glass Mylar, the new process becomes significant only 
at a sufficiently high electric bias field. In our experi- 
ments, we measured the dielectric constant of a 15 /jm 
thick Mylar film and a 300 /j,m thick BK7 wafer, both 
contained in a 3 He immersion cell. All measurements 
were performed in the linear response region and at the 
frequency / = 10 kHz. 

Figs. □ and show the excess dielectric response of 
Mylar after the temporary application of an electric bias 
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FIG. 1: Dielectric response of Mylar to a temporary applied 
bias field with maximal value F for a waiting time t w = 30 s 
and different sweep times at a temperature T = 10 mK. Solid 
lines are model curves explained in the text. The schematic 
illustrates the bias field sweep and the definition of t = 0. 

field for the waiting time t w — 30 s and at temperatures 
of T = 10 mK and T — 20 mK. The sketch in Fig. □illus- 
trates the evolution of the bias field sweep. Note that the 
zero point of time is defined to be at the completion of the 
sweep and, here, the time axes are logarithmic. In Fig.^ 
the sweep times t s vary between 0.3 s < t s < 300 s at elec- 
tric bias fields of maximal strength F = 6.7 MV/m, while 
in Fig. t s — 10 s for all curves at different maximum 
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FIG. 2: Dielectric response of Mylar as in Fig.0 but at con- 
stant sweep times for several maximal bias fields F. The data 
are taken at T = 10 mK (main figure) and T = 20 mK for the 
inset. Solid lines are model curves as explained in the text. 

bias fields 0.3MV/m < F < 8.7 MV/m. 

A temporary applied bias field causes relaxational pro- 
cesses that bring TSs out of equilibrium. In Q, strong 
deviations from the expectations due to relaxation were 
observed in the dependence of the waiting time (Fig. Q] 
in 0). These results were qualitatively explained by 
considering the evolution of energy splittings E(t) = 
v/A^ + (A + P F(i)) 2 while sweeping TSs through zero 
asymmetry A + pF(t) = (assembling an avoided level 
crossing situation) . If the combination of electric ac- and 
bias-field changes fast, systems tunneling slowly cannot 
adjust to their momentary equilibrium while being sym- 
metric, i.e. A + pF(i) < Ao. Such TSs are driven non- 
adiabatically out of equilibrium and, if part of a strongly 
coupled pair, suddenly contribute to the dielectric re- 
sponse. Due to the range of relaxation times, the ex- 
cess response decays in time as log(ro/t), as long as the 
ac-field dominates the driving effect, i.e. for small sweep 
rates F/t s . At faster sweep rates the bias-field drives TSs 
non-adiabatically. This situation causes deviations from 
logarithmic decays 0, E3- In Figs. UJ and the transi- 
tion between both non-adiabatic limits can be seen |ld| . 
Here, relaxational processes can be omitted, as compari- 
son to the waiting time dependence reported in shows. 
The slopes of the solid lines in both figures, representing 
logarithmic decays as expected for non-adiabatic ac-field 
driving, are calculated from the results in 0. They de- 
pend on temperature and maximal bias field but not on 
sweep rate and waiting time 0. The decay times tq, de- 
fined as the times when the model curves extrapolate to 
de/e — 0, depend on the many parameters of the model, 
but we take them here as fit parameters. 
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FIG. 3: Decay times for Mylar obtained from model curves 
of the previous Figures due to a model including dominant 
pair breaking by non-adiabatic driving of TSs l] . The solid 
lines are model curves without (and the dashed lines with) an 
additional contribution to the dielectric response that decays 
independent of temperature and bias field history. 



Fig. [31 shows double logarithmic plots of decay times 
as a function of the sweep rate F/t s . Solid points repre- 
sent the decay times extracted from measurements with 
a maximal bias field of F = 6.7 MV/m and different 
sweep times, for T = 10 mK (Fig. ^| and T — 20 mK in 
the upper and lower plot, respectively. Non-adiabatic 
driving causes tq to be a function of the sweep rate 
F/t s [l|,llO(. If bias field driving is dominant (fast sweeps) 
the model introduced in Q] proposes a linear depen- 
dence To = £,(T)t s /F. For dominant ac-field driving 
(slow sweeps) the exact form of the functional depen- 
dence is unknown. All data sets presented in Fig. [21 con- 
tain the transition between these two limits (fast and 
slow sweeps). The decay times extracted from our sweep 
time dependent data (solid points) show the same power 
law behavior over 4 orders of magnitude in F/t s , clearly 
indicating equal functional behavior in both limits. How- 
ever, we observe the relation To oc (t s /F) aA , in striking 
contrast to the expected linear dependence. 

The open circles in Fig. [21 present the decay times 
obtained from the bias field dependent measurements 
shown in Fig. [3 For bias fields greater than a crossover 
field F c ~ 1.5 MV/m we observe identical behavior as 
for the sweep time dependent measurements, that were 
taken at F — 6.7 MV/m, much greater than F c . Strik- 
ingly, for F < F c the behavior suddenly changes to the 
expected linear relation To cx t s /F (solid lines in Fig. [2J. 
The temperature dependence mirrors this situation: For 
F < F c we find t (T = 20mK) /t (T = lOmK) ~ 1.8 
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(solid lines), in agreement with the observed tempera- 
ture dependence of the decay times in the limit of dom- 
inant relaxational processes 0. For F > F c we ob- 
serve t (T = 20mK)/r (T = lOmK) ~ 1.8 04 (dashed 
lines), with the identical exponent as for the sweep rate 
dependence in this limit. Both the non-linear relation 
t oc (t s /F) 0A for F > F c as well as the temperature 
dependence cannot be explained by our model so far. 

Our BK7-sample is twenty times thicker than the My- 
lar film resulting in an accordingly weaker maximal bias 
field of F = 0.4 MV/m for our maximum applied volt- 
age. However, the electric coupling between TSs is about 
6 times stronger than in Mylar yielding similar excess 
responses in both samples for otherwise similar parame- 
ters. The sweep time dependence for BK7 is displayed in 
the inset of Fig.0]for t s < 10 s and in the main figure for 
t s > 10 s. Due to the stronger coupling between TSs in 
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FIG. 4: Sweep time dependence of the relative dielectric 
constant of BK7 after transient application of a bias field 
F = 0.4 MV/m at the temperature T = 20 mK for long 
sweep times ts > 10 s in the main figure and short sweep 
times t s < 10 s in the inset. 

BK7, for t s > t w a sizable decay of the excess dielectric 
response already happens while the bias field is still de- 
creased back to zero. In a simplifying approximation we 
assume an effective increase of the waiting time propor- 
tional to the sweep time. For the data shown in Fig. 0] 
the linear time transformation t — > t + 0At s produces a 
set of almost identical curves. This ability to compen- 
sate for the unwanted effect described above, is gratify- 
ing but not essential for our analysis of the decay times, 
which arc almost unaffected by the above time transfor- 
mation (due to the logarithmic time scale). The absence 
of any sweep time dependence in BK7 - in striking con- 
trast to the observations in Mylar - can be understood 
in terms of very different lower cutoffs of the density of 
states for TSs, -Emm, in the two samples. Only TSs with 



slow tunneling frequencies are driven non-adiabatically. 
Hence, we expect to observe sweep rate independent de- 
cays, if .E m in is larger than a crossover splitting, A 0c , be- 
low which non-adiabatic driving of TSs is pres ent. Tak- 
ing into account all of our measurements [lOj, we esti- 
mate E min /h > 1.6 mK for BK7 and E min /h < 40 fiK 
for Mylar, in agreement with conclusions drawn from 
a saturation of the dielectric constant at low tempera- 
tures For relaxational processes bringing pairs out 
of equilibrium, an upper limit for the decay time can be 
estimated. By taking only one phonon processes of TSs 
with E — k B T at zero field and E — pF > k B T at 
applied field into account and assuming a typical dipole 
moment of p ~ 1 D we find r < t w pF/k B T ~ 100 s. 
In contrast, we observe decay times scattered around 
To = 7100 ±500 s in our experiments (Fig.0J|. The excess 
dielectric response of BK7 is therefore primarily caused 
by a novel process, for which the decay time exhibits at 
most a very weak dependence on temperature |l0| . 

We will now show that the decay time dependence on 
the sweep rate To(F/t s ) in our Mylar sample (Fig. |3J) 
can as well be described by including such an addi- 
tional process with constant decay time. We assume 
that TSs contributing to the novel process do not con- 
tribute to non-adiabatic driving and therefore 5s /e cx 
[1 - a(F)} \n[£(T)t s /(Ft)] + a(F) ln[f /i]. Here £(T)t s /F 
and To are the decay times due to non-adiabatic driving, 
containing the temperature dependent pre-factor £(T), 
and the novel process, respectively, and a(F) is the am- 
plitude of the novel process. The decay times obtained 
from the logarithmic theory curves are then given by 

log (t ) = [1 - a(F)] log [£(T)t B /F] + a(F) log [f ] . (1) 

All lines in Fig. [3] are obtained from QJ, where we as- 
sumed a(F > F c ) — a Q and a(F < F c ) = 0, i.e. that 
the novel process disappears for F < F c , in accordance 
with experiments. The fit parameters are To ~ 4.7- 10 3 s, 
ao ~ 0.6 for both temperatures, £(10mK) ~ 620 MV/m 
and £(20mK)/£(10mK) ~ 1.8. By including the novel 
process we can now recalculate the model curves for 
Figs, ^and [3 as well as for the data in a previous pub- 
lication Hereby, the slopes of all model curves only 
depend on two commen fit parameters, the dipole mo- 
ment of p = 1.2±0.1 D and the interaction constant, 
PqUq ~ 6.3 • 10~ 4 , where Pq is the density of states of TSs 
and Uq the mean interaction energy times volume. Note, 
that due to the influence of the novel process, our value 
of PoUo, is slightly higher than previously reported 0. 
Using logarithmic decay curves for the novel process we 
find good agreement with our data by assuming the same 
dependence of the slopes of the decay curves on tem- 
perature and maximal bias field as already proposed in 
the "dipole gap" theory y| and observed in our previ- 
ous experiments for relaxational pair breaking and non- 
adiabatic driving processes 0. Thus, the slope of the 
decay curves (in their logarithmic region) is governed by 
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the general dependence on temperature of the dielectric 
response as well as how far the bias field drives the TSs 
out of equilibrium. In order to rea ch g ood agreement 
with all our measurements on Mylar [lOj, we necessarily 
must make the assumption in Q that the novel process 
excludes additional contributions to the excess dielectric 
response due to other processes, i.e. non-adiabatic driv- 
ing. Undoubtedly, if the dynamics of strongly coupled 
pairs of TSs governs the other contributions, the novel 
process influences clusters as well. 

The temperature independence of tq points to a pro- 
cess involving quantum mechanical tunneling, wherein 
thermal transitions are unimportant. The general as- 
sumption of a flat distribution of tunneling parameters 
p(A) = po in a disordered solid yields the distribution of 
tunneling times p(t\) = Po/t\, with t\ oc exp(A). As- 
suming that strongly coupled pairs stop contributing to 
the dielectric response once they are in thermal equilib- 
rium, as in the " dipole gap" theory, we find a logarithmic 
decay Se/e oc f^ Km " dTp(T)Q(t — r) = p hi(To/i), with 
= t and the Heavyside step function Q(t - r). A 
strong electric bias field couples to dipole moments, gen- 
erally present in amorphous solids due to displacements 
of atoms, and leads to structural rearrangements of the 
atoms changing the potential energy distribution. This 
situation yields a different set of low energy states than 
the TSs provided at zero field. Nevertheless, in thermal 
equilibrium, experimental observations are identical be- 
cause of the glassy distribution of parameters. Primarily, 
we are interested in the effect of a bias field sweep on TSs 
contributing at zero field to the dielectric constant. 

Simulations suggest that the additional low energy in- 
ternal degrees of freedom observed in glasses include co- 
herent movements of the order of 20-100 atoms [llj . The 
interaction of that many atoms yields a multiple level 
scheme, wherein a TS can contribute to the excess di- 
electric response, if the two lowest levels are close enough 
in energy, e.g. A < k^T and E ~ k^T. A strong bias 
field pF 3> k-gT couples to dipole moments and reor- 
ganizes the energy eigenstates present in the glass. It 
thereby brings TSs, including strongly coupled pairs, out 
of thermal equilibrium. Moreover, energy states that are 
thermally not accessible at zero bias field will be occupied 
during the field sweep by quantum mechanical tunneling. 
It is conceivable that successive tunneling processes while 
the bias field is swept tend to bring TSs into metastable 
states that are lower in energy than their neighbor po- 
tential wells (in configuration space) during the entire 
field sweep, i.e. tunneling processes tend to be unidirec- 
tional in a non-equilibrium situation. Once back at zero 
bias field, these TSs are generally still trapped in their 
metastable states and separated from the lowest energy 
states (at zero bias field) by a large potential barrier. In- 
termediate states temporarily occupied during the field 
sweep are now statistically higher in energy and ther- 
mally not accessible (SE ^> fcsT). While metastable 



states in single wells are occupied that cannot be counted 
as TSs, additional metastable two level states exist. Be- 
cause of the flat distribution of asymmetry energies, the 
overall number of TSs contributing is the same at any 
time. However, TSs out of equilibrium that belong to 
strongly coupled clusters contribute to the excess dielec- 
tric constant as being isolated. Our interpretation in 
terms of structural rearrangements should be considered 
as a rough outline of a possible scenario. It calls for a de- 
tailed examination, including the influence of metastable 
conditions on the dielectric response dependent on the 
broadly distributed time constants of the different in- 
volved processes. 

We observe the novel process in BK7 at relatively small 
bias fields. However, a minimum bias field is required 
in the polyester glass Mylar possibly due to a potential 
barrier necessary to rearrange hydrocarbon chains. 

In conclusion, after the temporary application of a suf- 
ficiently large electric bias field to a glass sample at ultra 
low temperatures, we find a long-lasting contribution to 
the excess dielectric response that decays logarithmically 
in time with a decay time independent of temperature 
and bias field history. We assign this novel process to 
field -induced structural rearrangements in the glass that 
decay slowly by quantum mechanical tunneling through 
large potential barriers. Strongly coupled pairs involved 
cause an excess dielectric response as long as they are not 
in thermal equilibrium. 
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